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Abstract 

Braneworld theory provides a natural setting to treat, at a classical level, the cosmological effects of vacuum 
energy. Non-static extra dimensions can generally lead to a variable vacuum energy, which in turn may explain the 
present accelerated cosmic expansion. We concentrate our attention in models where the vacuum energy decreases 
as an inverse power law of the scale factor. These models agree with the observed accelerating universe, while 
fitting simultaneously the observational data for the density and deceleration parameter. The redshift at which 
the vacuum energy can start to dominate depends on the mass density of ordinary matter. For fl m — 0.3, the 
transition from decelerated to accelerated cosmic expansion occurs at zt ~ 0.48 ± 0.20, which is compatible with 
SNe data. We set a lower bound on the deceleration parameter today, namely q > — l + 3f2 m /2, i.e., q > —0.55 for 
= 0.3. The future evolution of the universe crucially depends on the time when vacuum starts to dominate 
over ordinary matter. If it dominates only recently, at an epoch z < 0.64, then the universe is accelerating today 
and will continue that way forever. If vacuum dominates earlier, at z > 0.64, then the deceleration comes back 
and the universe recollapses at some point in the distant future. In the first case, quintessence and Cardassian 
expansion can be formally interpreted as the low energy limit of our model, although they are entirely different 
in philosophy. In the second case there is no correspondence between these models and ours. 
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1 Introduction 



It is now widely accepted that the present universe is accelerating and spatially flat. Evidence in favor of this is 
provided by observations of high-redshift supernovae la pDh|I3, as we ^ as other observations, including the cosmic 
microwave background and galaxy power spectra |7]-|13|. 

Since ordinary matter is gravitationally attractive, one can conclude that the source of cosmic acceleration is some 
kind of unknown matter, which violates the strong energy condition. According to dynamical mass measurements, 
the total amount of ordinary matter in the universe, including dark matter, can only account for 30% of the critical 
density. Thus, the other 70% corresponds to a mysterious form of matter, usually called dark energy, which remains 
unclustered on all scales where gravitational clustering of ordinary matter is seen. 

Therefore, in the past few years an active field of research has been the discovery of models of the universe 
in which the expansion is accelerating. The simplest candidate for dark energy is the cosmological constant ku\ 
|14j-|15j. In this approach, A(4) is introduced by "hand" as a parameter in Einstein's theory of gravity. However, if 
A( 4 ) remains constant one faces the problem of fine-tuning or "cosmic coincidence problem" |lfi| . which refers to the 
coincidence that p vac and p m are of the same order of magnitude today. 

A phenomenological solution to this problem is to consider a time dependent cosmological term, or an evolving 
scalar field known as quintessence JU-^H]. A variable Am, as well as quintessence, can be modeled as the energy of 
a slowly evolving cosmic scalar field <fi with an appropriate self-interaction potential V(4>) to account for the evolution 
of the universe. The abundance of quintessence models is due to the fact that for any given scale factor a(t) and 
some known forms of energy density Pknown(t) (made of radiation, matter, etc), it is always possible to find a V((j>) 
that explains the observations [2*0"| . 

Alternative explanations for the acceleration of the universe, beyond dark energy, include phantom energy |21j- 
|24|. certain modifications of general relativity US]- [33], the gravitational leakage into extra dimensions |34|-|36|. 
Chaplygin gas dmUS] as well as Cardassian models @H]-@I]. More recently, in a series of papers Vishwakarma [4*2"] . 
[4"3"| and Vishwakarma and Singh 02] argue that present observations can successfully be explained by a decelerating 
model of the universe in the mainstream cosmology, without invoking dark matter or dark energy. 

In the present work, we study the accelerated cosmic expansion as a consequence of embedding our observable 
universe as a brane in extra dimensions. Our motivation is that braneworld theory provides a natural setting to 
treat, at a classical level, the cosmological effects of vacuum energy. Firstly, the theory links the vacuum energy to 
the fundamental quantities A(4) and G. Secondly, in this theory the vacuum energy is completely determined by the 
geometry in 5D through Israel's boundary conditions Said another way, in brane theory the cosmological 

term is determined by the bulk geometry and cannot be put by hand as in Einstein's theory. Thirdly, non-static 
extra dimensions can generally lead to a variable vacuum energy 0E|, and consequently variable A/ 4 ) and G, which 
in turn may explain the present accelerated cosmic expansion. 

In fact, in recent papers |49 | -|50 | we showed that braneworld models with variable vacuum energy (and zero 
or negative bulk cosmological constant) agree with the observed accelerating universe, while fitting simultaneously 
the observational data for the density and deceleration parameter. In those papers we were mostly interested in 
the general behavior of the models. We obtained precise constrains on the cosmological parameters as well as 
demonstrated that the "effective" equation of state of the universe can, in principle, be determined by measurements 
of the deceleration parameter alone. 

In this paper we are interested in more subtle details. For example, we would like to know the rcdshift of transition 
from deceleration to acceleration. Namely, based on our model, can we predict, or at least narrow down, the redshift 
interval for the transition? What is the time at which the vacuum energy should start to dominate in order to explain 
the observed cosmic acceleration? How is it related to the current mass density of the universe? Does the present 
deceleration parameter depend on it? What is the future evolution of the universe? Will it continue to accelerate 
forever? 

In the scenario discussed here, the acceleration of the universe is related not only to the variation of vacuum 
energy and cosmological term, but also to the time evolution of G and, possibly, to the variation of other fundamental 
"constants" as well. 

The paper is organized as follows. In section 2 we give a brief summary of the equations for homogeneous 
cosmologies in 5D based on braneworld theory. In section 3 we show how to incorporate a varying vacuum energy 
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into the scheme. Observational constraints on the model are discussed in section 4. In section 5 we discuss the 
transition from deceleration to acceleration. 

We show that the crossover point strongly depends on the energy density of ordinary matter. The possible values 
of the redshift of transition are spread over a finite interval, regardless of the moment in time when vacuum starts 
to dominate over ordinary matter. This is a clear indication that the effects of vacuum rapidly decrease with the 
increase of z. We also show that the value of the deceleration parameter today is bounded bellow. 

We discuss the cosmological parameters for the case where the matter density today is one-third of the critical 
density. We find that the future evolution of the universe crucially depends on the time when vacuum starts to 
dominate over ordinary matter. Namely, if the vacuum contribution starts to dominate only recently, at an epoch 
z < 0.64, then the universe is accelerating today and will continue that way forever. But, if the vacuum dominates 
earlier, at z > 0.64, then the deceleration comes back and the universe recollapses at some point in the distant future. 

Finally, in section 6 we give a summary and discussion. We show that the original Cardassian expansion, proposed 
by Freese and Lewis |40) . as well as dark energy models (with constant wx) can be interpreted as the low energy 
limit of our ever-expanding models. 



2 Homogeneous cosmology in 5D 

In order to facilitate the discussion, and set the notation, we start with a brief summary of the pertinent ideas and 
equations in the braneworld scenario. In this scenario our homogeneous and isotropic universe is envisioned as a 
singular hypersurface embedded in a five-dimensional manifold with metric 



dS 2 = n 2 (t, y)dt 2 - a 2 (t, y) 



dr 2 



r 2 {d6 2 



* 2 (t,y)dy 2 , 



(1 - kr 2 ) 

where t, r, 9 and <f> are the usual coordinates for a spacetime with spherically symmetric spatial sections and k 
0, +1, — 1. The metric is a solution of the five-dimensional Einstein equations 



(1) 



(5) G 



AB 



^RAB-l 9 AB ( - 5) R=kf 5) ^T AB , 



(2) 



where ^'Tab is the five-dimensional energy-momentum tensor and k^ is a constant introduced for dimensional 
considerations. 

The energy- momentum tensor on the brane t^ v is separated in two parts, 



(3) 



where a is the tension of the brane in 51?, which is interpreted as the vacuum energy of the braneworld, and T^ v 
represents the energy-momentum tensor of ordinary matter in AD. 

There are two assumptions relating the physics in AD to the geometry of the bulk. They are: (i) the bulk 
spacetime possesses Z2 symmetry about the brane, and (ii) the brane is embedded in an Anti-de Sitter bulk, i.e., 
^'Tab is taken as 

^T AB = A {5) g AB , (4) 

where A( 5 ) < 0. 

As a consequence of the first assumption, the matter and vacuum energy density in AD become completely 
determined by the geometry in 51? through Israel's boundary conditions. Namely, for a perfect fluid 



= (p + p)u^u v - pg^, 



with 

the matter density is given by [48) 



P = IP, 



(-2) 



*?6)(7 



a 
a 



11 

n 



(5) 
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(7) 
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and the vacuum density is 

(-2) 



a! ri 



(37 + 2)- + - 

a n 



(8) 



brane 

where a prime denotes derivative with respect to the extra variable y. 

The second assumption has two important consequences. Firstly, since ^T^i = 0, it follows that the energy- 
momentum tensor on the brane t^ u is a conserved quantity, viz., 



r 



= 0. (9) 



Secondly, the field equations J2J admit a first integral, namely, 

fc (5) A 0) , (a'\ 2 k C 



. 2 1,2 

a 



an ) 6 V cl& 



-,2 1 „4- 



(10) 



where C is a constant of integration which arises from the projection of the Weyl curvature tensor of the bulk on 
the brane. Evaluating (|10|l at the brane, which is fixed at some y — ybrane — const, as well as using Q and (jSJ), we 
obtain the generalized Friedmann equation, viz., 



k%\ 9 3fc 3C 



3( ^j =A (4)+ 8^ + ^y-^, 



where a (t) = a(t,y bra ne), and 



1 ( k 2 a 2 \ 
A(4) = 2 fc (5)( A (5) + ^-)^ (12) 

kt-,0- 

8ttG=^-. (13) 

These quantities are interpreted as the net cosmological term and gravitational coupling in 4 dimensions, respectively. 

Equation Ijlll) contains two novel features; it relates the fundamental quantities A( 4 ) and G to the vacuum energy, 
and carries higher-dimensional modifications to the Friedmann-Robertson- Walker (FRW) cosmological models of 
general relativity. Namely, local quadratic corrections via p 2 , and the nonlocal corrections from the free gravitational 
field in the bulk, transmitted by the dark radiation term C/a 4 . 

Except for the condition that n = 1 at the brane, the generalized Friedmann equation is valid for arbitrary 
Q(t, y) and n(t, y) in the bulk 07]. This equation allows us to examine the evolution of the brane without using any 
particular solution of the five-dimensional field equations. In what follows we will omit the subscript 0. 



3 Variable vacuum energy 

In equation G and A( 4 ) are usually assumed to be "truly" constants. However, the vacuum energy density a 
does not have to be a constant. From JHJ) it follows that a depends on the details of the model. Indeed, we have 
recently shown |48| that there are several models, with reasonable physical properties, for which a variable induces 
a variation in the vacuum energy a. 

3.1 Variable vacuum: an example 

As an illustration, let us consider the class of solutions to the field equations J2Jl for which the metric coefficients in 
Q are separable functions of their arguments. For this class of solutions, without loss of generality we can set 

n = n(y), a(t,y) = a(t)Y(y), $ = $(*). (14) 



4 



From G04 = it follows that 

TrJ-<(y). 1^-4 <»> 

where £ is a separation constant. Consequently, 

*(t) = Aaf 1 -^, (16) 

where A is a constant of integration. Thus, for any £ ^= 1, $ is a variable function of i. In what follows, for simplicity 
of the notation, instead of ( we will use (3 = — (£+2)/3 (or ( = -3/3-2). With this notation we have $(f) = Ad 3< - j3+1 \ 
Substituting this expression into © we find 

- D <"-® ( i7) 



where we have introduced the constant D = +l)Y^ rane Y 3 fJ^ / A(j +l)k^ 5 y This epitomizes the general situation 
where a is a function of time. It is worth noting that only the assumption of separability of the bulk metric (fT|) 
underlies equation (|17|) . i.e., we do not need to know the details of the solution in the five-dimensional bulk. 

From a physical point of view, the vacuum energy (|17fl implies (G/G) = —3(j3+l)H. Conversely, if we extrapolate 
the present limit |G/G| = \g\H, and assume that g is a constant 1 , then we obtain a cosmological model where the 
vacuum energy is given exactly by (|17(l . with g = — 3(/3 + 1). This model is consistent with the observed acceleration 
and flatness of our universe • 

3.2 Our model 

In general, for variable vacuum energy, the conservation equations for a perfect fluid which satisfies ®, yield 

p+3p(j + l)-=-&, (18) 
a 



From which it follows that 

C 1 



a 3(7+l) a 3 (T +1 ) 



a 3 ^ +1 W, (19) 



where G is a constant of integration. For the case of constant a, we recover the familiar relationship between the 
matter energy density and the expansion factor a, viz., 

The second term in l|19|l is the contribution associated with the variation of vacuum. The variation of the vacuum 
energy is deeply rooted in fundamental physics. The simplest microphysical model for a variable <r, as well as for 
A( 4 ) and quintessence, is the energy associated with a slowly evolving cosmic scalar field <f> with some self-interaction 
potential V(<p) minimally coupled to gravity. The potentials are suggested by particle physics, but in principle V{<f) 
can be determined from supernova observations |51 j -|53 | . 

In this paper, instead of constructing a field theory for the time evolution of the vacuum energy 2 , we employ our 
previous example as a guide. Namely, if during the expansion of the universe a decays as in (|17f) . then from 119|) it 
follows that 

- - - 

P ~ ^3(7+17 + a 3(/3+l) ' ( > 

1 The physical meaning of this assumption is that the variation of g is much "slower" than that of H and G. Namely, \g/g\ << \H/H\, 
\g/g\«\G/G\ 

2 According to Padmanabhan 20; it is trivial to choose the "appropriate" potential V(<j>) such that we can explain the observations, 
for any given pair a(t) and p(t) 
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Now, from the conservation equation (|18|) we get 



gfr - (3) 

,J = fTo+ (/3+l) a 3(/ 3+ i)' ( 22 ) 

where (To is a constant of integration. We will assume this form, with ^ 0, for the vacuum energy 3 . It reduces to 
the case of constant vacuum energy a = ero, f° r P = 7- However, for [3 ^ 7 it generates models with variable er, G 
and A(4). We immediately notice that the positivity of G — kf 5 \<r/6 requires (3 < 7. 

3.2.1 The effective density 

The total energy density of the universe, ptotai = a + P, can be written as 

Ptotal = CO + p e //, (23) 

where we have introduced the "effective" density p e ff as 

= ^+1) + [j^J ) ^TT) • ( 24 ) 

This effective density is the one that drives the evolution of the universe. Indeed, in the present model the generalized 
Friedmann equation becomes 



2 1 / 1,2 \ ui _„ 1.4 




k (5) a0 , fc (5) 2 3fc 3C 



3 - ^'i "1 ( -V.-., ^ 1 + ^-/'. n 1 T7''' ' ' " - - 

We note that, distinct from l|21|) . the vacuum contribution in the effective density ()24|l is multiplied by the factor 
(7 + l)/(f3 + I), which is larger than 1, because (3 < 7. 

In what follows we will set 7 = 0, in view of the fact that our universe is matter-dominated (p = 0). In addition, 
we can set k = 0, because astrophysical data from BOOMERANG [HI] and WMAP indicate that our universe is 
flat. Also we can consider C = 0, since the constant C, which is an effective radiation term related to the bulk Weyl 
tensor, is constrained to be small enough at the time of nucleosynthesis and it should be negligible today. Also, in 
order to avoid an exponential expansion of the universe in its asymptotic limit, we assume 

k 2 

A (5) + "f^o = 0. (26) 
Thus, in the case under consideration the evolution of the universe will be governed by 

a) -—Q-P'Jt + l2 P °"> (2?) 
with 

_ C D/W+l) 

Mf ~ ^ + a 3(/3+D ■ (28) 

The cosmological term A( 4 \ is not constant, but evolves according to 

(4) 6 (/3+I)a 3 ^ +1 ) 12 + { V > 



'With this assumption g is not constant. 
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3.2.2 Asymptotic behavior 



For (3 = 7 = 0, we recover the usual picture, i.e., p = p(l + z) 3 , a = o~q and A( 4 ) = 0. For (3 < 7 = 0, the vacuum term 
is initially negligible, which means that p e // approaches the typical matter density in FRW models and (a / 1 p) — ► 0. 

If Co is positive, then the universe is in continuous expansion. When the vacuum term in i|28[) is so large that the 
ordinary matter contribution can be neglected, we find 



a 



3(/3+l) ^ (5) u yh 1 .1 



t 2 . (30) 



The corresponding "deceleration" parameter q = —aa/a is given by 



g~ 2 . (31) 

which indicates that the expansion is accelerated for {3 < —1/3. At this stage of the evolution G is constant and the 
cosmological term A/4) varies with time. Namely, 

3.2.3 The vacuum energy takes over 

We note that, for j3 < 0, the first term in (|28[) decreases in time faster than the second one. Therefore, D has to be 
chosen in such a way as to make the second term in (|28|) important at the right time to explain the observations. 
In order to do this, we find useful to introduce the auxiliary quantity z eq . This is the redshift at which the vacuum 
and matter components in (|28J) become equal to each other 4 , i.e., C/a 3 = D/{[3 + l)a 3 ( /3+1 ) at a = a(z eq ). Since 
a = a/(l + z), where a is the present value of a, we obtain 

c =4w- (33) 



The appropriate D follows from the evaluation of p e f f today. We find 

D (fi + l)p eff 



a 3 ^ 1 ) [! + (! + z eq )W] ■ 



(34) 



where p e f t is the present value of the effective density. 

We note that the case where D = and C 7^ is formally obtained from our equations by setting z eq — — 1 , for 
any (3 < 0. Similarly, the case where C = and D ^ is formally attained in the limit z eq — > 00. 



4 Observational constraints on (3 



Although the evolution equations (|27|I - H28|I contain four constants: C, D, oo, and fc 4 5 ^, there are only two parameters 
in the model, viz., (3 and z eq . The aim of this section is to find out the physical restrictions on the parameter (3. 

With this goal, we express the relevant quantities in terms of these parameters and the present-value of the 
density of ordinary matter p m — C/a 3 . Let us start with the effective density H28JI . Using l)33|l and Q34J1 ■ it can be 
written as 



Peff = Peff 



(l + z eq ) 



[! + (! + z eq ) 



3/3 



il + zf 



1 



1 + z 
1 + z eq 



3/3' 



(35) 



4 In order to avoid misunderstanding: the parameter z eq 
denote as zt- 



is not the redshift of transition from deceleration to acceleration, which we 
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The ratio of effective density to ordinary matter density is given by 

Peff _ 1 | ( 1 + z x '" 

Pm \ 1 ~h ^eg 



(36) 



Thus, at the present time 

Peff = p m F(z eq ,(3), (37) 

where 

[1 + (1+^] 

F(z ^' /3) - (TT^F ■ ( } 

With this notation, we have C = a 3 p m , D — a 3 ^ +1 ^(F — l)(/3 + l)p m . Consequently, 

Pe ff -p m [l + (F-l)(l + z) 3/3 ]. (39) 

Thus, in our formulae the case of constant vacuum energy, for which D = and /0 e // = p m , corresponds to F = 1, 
for any /3 < 0. For large redshifts p e // ~ Pm, while at the present time peff = PmF. We note that F can be very 
large for large values of the parameter Z eq ■ 

4.1 Positivity of G 

We now proceed to calculate oo- Evaluating (|27[) today, and using l|37|) . we have 

1 



kf 



3H* ^ 



a p m F(z eq ,/3) + -p m F {z eq ,0) 
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The constant fcfu is given by fcfu = 487rG/er. Thus, using i|22[l . with 7 = 0, l|34|l and (|38|l . we obtain 



(40) 



4 48.G 

fc(5) " [a -/?p m (F-l)]' (41) 
where G is the present value of the gravitational "constant" G. Feeding this expression back into H40(l we obtain 

p m F*/2 + f}(F-l)] 

^o^l ' (42) 

where fi m is the present value of the mater density parameter fi m = 8irGp m /3H 2 . The vacuum term (1221) can be 
written as 

a{z,Ze q ,l3)=a -l3p m (F-l)(l + z) 3 ^ +1 K (43) 
Evaluating this equation today, we have 

-i ro \ - n m FlF/2 + (3(F-l)} 

a(z eq ,P,il m ) = p m - (44) 

1 — \l m r 

We now calculate the constant kf 5 y For this we substitute l|42|l into (|41|l . We get 

k (5) {z eg ,P)- _^ [F/2 + p{F _ 1)]F - ( 45 ) 

We are now able to express general physical conditions on (3. Since 8nG — fcj^er/6, the positivity of G demands 
kf 5) a > 0. Thus, for any given f2 m and z eq , the allowed values of j3, in (|44|l and (|45|) . are those for which 5 

^ 5) (z e «,/3,^m) > 0, a(z eq ,(3,Q m ) > 0. (46) 



5 We note that a > <r. Also, the possibility k^(z eq , f3) < 0, and a(z eq , (3,Q m ) < is excluded by the fact that fc? g . in i s a rea l 
number. 
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4.2 G/G 

More stringent restrictions on j3 follow from observational constraints on the variation of G. In terms of the Hubble 
parameter, the time evolution of G is usually written as (G/G) — gH, where g is a dimensionless parameter. In our 
model we have 

<7 = -3(/3+l)(l-^). (47) 

We note that 3 = for (3 = 7 (i.e., (3 = for ordinary matter). From (|42(l and (|44Jl . we obtain the present value of g 
as 

g(z eg , (3, n m ) n m F[F/2 + 8(F-l)] (48) 

The abundance of various elements as well as nucleosynthesis are used to put constraints on g today. The present 
observational upper bound is 6 

\g\ < 0.1. (49) 

Thus, the choice of (3 has to guarantee the fulfillment of this condition. It clearly pushes the values of (3 either toward 
/3 « or /? -1. 

4.3 Deceleration parameter q 

In order to consider other observational constraints on /3, let us introduce the deceleration parameter 

, = _g = _ 1+ <™ M 

which in the case under study becomes 

q{) ~ 1+ (2,0+Mf) [1 + (F-1)(1 + ^] ' (51) 
where p e f / (z) and <7o are given by H36|l and (|42|l , respectively. Evaluating this equation today, we have 

q[z eq ,p,U m )- 1 + d F [F + 2B(F-1)] ' (& j 

In the case where (3 = 7, the present-day acceleration l|52|l reduces to 

q = 2-Ztt m /2. (53) 

In particular, for Q, m = 1 we obtain q = 1/2 as in the dust FRW cosmologies. We note that (|53|l is positive for 
any physical value of fi m , which means that a brane- universe with constant vacuum energy must be slowing down 
its expansion 7 . However, for j3 ^ 7, this is no longer so; the vacuum energy a is now a dynamical quantity which 
changes this picture. 

The choice of (3 in l|52|l has to be consistent with recent measurements which indicate that the current universe 
is speeding up its expansion with an acceleration parameter which is roughly 

q = -0.5 ±0.2. (54) 



6 A comprehensive and updated discussion of the various experimental and observational constraints on the value of g (as well as on 
the variation of other fundamental "constants" of nature) has recently been provided by Uzan 1551 

7 This contradicts the observational fact that the universe is speeding up and not slowing down 1541 . 
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4.4 The cosmological constant 

Finally, for the cosmological parameter Am at the present day we have 

A(4) =Z(ze q ,f]^ m )H 2 1 (55) 

where 



t, ao \ 3/3(f - l)[q»F 2 + P(F - 1)(1 + n m F)} 

e(^,/3,0 m ) = F [F + 2f3(F-l)]) ' (56) 



Thus, A(4j = for (3 = 7, as expected. Otherwise, for (3 7^ 7 = 0, Ku\ has to be positive in order to explain the 
present acceleration. 



5 The auxiliary parameter z eq 

We now turn our attention to parameter z eq . Firstly, we confine the range of zt, the redshift of transition from 
deceleration to acceleration, which is a solution of q(zx) — 0. Secondly, we set a lower bound on the value of the 
deceleration parameter today. Thirdly, we study in some detail the models with Cl m = 0.3. An interesting discovery 
here is the possibility of recollapsing models. In what follows we will only select (3 from the range of values allowed 
by the requirements jS| , jH^I and (|51Jl . 



5.1 Transition from deceleration to acceleration 

There is plenty of observational evidence for a decelerated universe in the recent past, see e.g. [2>]-[57|. However, 
the dominance of the vacuum energy at some z eq > does not guarantee the present acceleration of the universe. 
For this, the vacuum energy has to dominate long enough as to overcome the gravitational attraction produced by 
ordinary matter. However, for every fixed value of f2 m , we can find a z™ m for which the transition from deceleration 
to acceleration occurs at some zt (0 < zt < z™" 1 ). 

The transition is then guaranteed for z eq > z™" 1 . A straightforward numerical study of (|51|l . in the interval 
(z™",oo), reveals that zt is bounded above. Thus in our model, for every Q m the transition from deceleration to 
acceleration occurs in a finite redshift interval 

z^ ln <z T < z T lax , (57) 

where the precise value of the lower and upper bounds depends on the density parameter for ordinary matter. The 
existence of an upper bound is an evidence that the dynamical influence of vacuum energy rapidly decreases for 
redshifts z > z T lax . We now proceed to show our numerical results for various values 8 of £l m . 



&m = 0.1: We find that a transition from deceleration to acceleration is only possible if the vacuum energy starts 
to dominate in an epoch before z eq > z™ q n — 1.3, otherwise the universe would be still today in a deceleration phase. 
The actual redshift of transition is zt = 1.1. A detailed investigation of q(zx) = 0, with q given by (|51|l . shows that 
zt < z T lax = 1.81 for z eq in the range 1.3 < z eq < 00, i.e., the transition occurs in the interval zt = (1-10, 1.81). Thus, 
for the whole range of z eq , we find that the redshift of transition from deceleration to acceleration, for Cl m = 0.1, is 9 

zt ps 1.46 ± 0.36. (58) 

8 A reliable and definitive determination of O m has thus far eluded cosmologists. However, the observational data indicate that 
fip ~ 0.1 — 0.3 seem to be the most probably options. 

9 We note that these numbers are approximate; they depend on the specific choice of in the range that satisfies the conditions 141 i I 
and 1491 . However, since \Af3//3\ is small and decreases with the increase of z eq (for an illustration see Table 1), the variation of these 
numbers is negligible and does not change the picture here. 
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= 0.2: The transition is guaranteed for z eq > 0.7. It occurs in the redshift interval zt = (0.58,1.10), or 
equivalently 

z T « 0.84 ±0.26, (59) 

and any value of z eq > 0.7. 

= 0.3: The vacuum contribution must start dominating at an epoch earlier than z eq = 0.4, otherwise there 
would be no enough time for a transition from a decelerating phase to an accelerating one today. For any z eq > 0.4, 
we find that the transition occurs at 

z T » 0.48 ± 0.20. (60) 

This is consistent with the value zt ~ 0.5 given by Turner and Riess |S7], and a little less than zt ~ 0.73 provided 
by Perlmutter 10 et al 0. 

The above discussion illustrates three things. Firstly, the fact that the redshift z eq , at which the vacuum energy 
starts to dominate, depends on the mass density of ordinary matter; in a low matter density universe the vacuum 
energy starts to dominate before than in a universe with high matter density. Secondly, that the redshift of transition 
increases with the decrease of the matter density (|58[1 - H60[) . Thirdly, that the vacuum effects only become important 
at the present epoch, making the transition a recent phenomenon. 

5.2 Lower bound on present deceleration (Upper bound on acceleration) 

The deceleration parameter today, which is given by 1(52(1 . is an increasing function of f3, for any fixed values of z eq 
and Q m . This means that its minimum value q m in is attained in the limit — ► — 1. 

Now, fixing Q m and selecting an appropriate (3 (as to satisfy the conditions discussed in the previous section) 
we find that q m in decreases for large values of z eq . As an example consider the values presented in Table 1; q m in 
increases for z eq < 0.64 and decreases for z eq > 0.64. Therefore, the lowest value of q m in is attained in the limit 
z eq — ► oo. Thus, from 1(52(1 we get 

g>-l + |o m , (61) 

for any value of z eq . For Cl m = 0.3, we find q > —0.55, which is consistent with recent observations 1)54(1. 

A similar analysis of (|55|l reveals that A(4), the cosmological "constant" today, is bounded above. Namely, 

A (4) < 3tf 2 (l - n m ), (62) 

for any z eq . Thus, if we take Q, m = 0.3, we get Am) < 2.1H 2 . However, this is just a statement that in the present 
model £l m + ^a (4) < 1 in view of the quadratic correction in the generalized Friedmann equation. 

5.3 The model for Cl m = 0.3 

Current dynamical mass measurements suggest that the matter content of the universe adds up to 30% of the critical 
density 11 . According to (|60|l the transition from deceleration to acceleration occurs at a redshift 0.28 < zt < 0.68, 
which confirms the idea that the accelerated expansion of the universe is a recent phenomenon. 

Thus, in our model, stars and galaxies with redshifts larger than z « 0.68 should reflect the kinematics of a 
decelerating expansion. This is compatible with galaxy formation, which can only take place if the gravitational 
attraction dominates a sufficiently long epoch over vacuum repulsion 12 . It also fits the observations of SN 1997ff at 
z « 1.7. This is the oldest and most distant type la supernova (SN la) discovered so far and provides direct evidence 
that the universe was decelerating before it began speeding up 60 . Other high redshift SNe known at z = 1.2 (e.g., 

10 It is encouraging that completely different models provide similar values for the redshift of transition from deceleration to acceleration. 
See for example 1581 

11 Radiation 0.005%, ordinary luminous baryonic matter 0.5%, ordinary non-luminous baryonic matter 3.5% and exotic (non-baryonic) 
dark matter "observed" through their gravitational effects 26%. 

12 For an updated analysis of the influence of dark energy on the first epoch of galaxy formation see, e.g., 1591 
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SN 1999fv, SN 1998eq), z = 1 (e.g., SN 1997ck, SN 1999fk) and z = 0.83 (e.g., SN 1996cl) can provide a direct test 
for deceleration at the time of their explosion. 

Also, from i|60|) it follows that galaxies with redshifts less that z w 0.28 should show evidence of an accelerating 
universe. Very-low redshift supernovae are crucial for reducing the uncertainty of the contemporary expansion rate. 

We now proceed to study in more detail the parameters of the model for Q TO = 0.3, which is favored by observa- 
tions. First of all, let us simplify the notation. To this end, we introduce the dimensionless parameters r]k and r) ao 
as 

K (5)Pm _ Op 

ism' Vao = 



(63) 



in terms of which the generalized Friedmann equation (|27|l becomes 



H 2 = H 2 Vk 



Peff 



Peff 



(64) 



We note that 



VkVa 



Vk 



F 



(65) 



which follows from (|42l) and 1451) . so that H = H today, as expected. 

Examination of 1441) . (|45|l and (|48|l reveals that, for every given z eq and Q m , the adequate values of (3, that satisfy 
the conditions 146|) and l|49l) . are spread over a small range. In that range q, rjk and r\ (7Q increase with f3, while a 
and A( 4 ) decrease with (3. For example, if we take z eq = 1 and fi m = 0.3, we find 13 (3 = —0.989 ± 0.011. The 
corresponding q increases from q = —0.46 for « — 1 to q = —0.35 for [3 — —0.978. The other quantities undergo a 
relatively smaller change, viz., 



m - 056+ 002 n - -2 3Q0+ 049 — - 5 343+0.196 . _ -, q2 Q+o.087 

Pm tl A 



(66) 



For the redshift of transition from deceleration to acceleration we find 



zt — (0.335, 0.298, 0.279) for (3 = (-0.999, -0.985, -0.978) respectively. (67) 

In Table 1, we illustrate the cosmological parameters for various values of z eq and Cl m = 0.3. For the sake of 
simplicity, we omit their small change over the range of (3 and present their mean values only , corresponding to 
the average of (3. Like we have said before the vacuum contribution must start dominating at an epoch earlier than 
z eq = 0.4, otherwise the universe would be still in a decelerating phase. 



Table 1: Cosmological parameters for Q, m — 0.3 fitting \g\ < 0.1 





-1< /3 < 


Vk 


(v/pm) 


q 




Ht 


(A (4 )/ff 2 ) 


z T = z q=0 ) 


0.4 


-0.840 


0.030 


10.052 


(-0.49,-0.38) 


7.724 


0.743 


1.853 


0.378 


0.5 


-0.949 


0.063 


4.791 


(-0.46,-0.35) 


1.602 


0.691 


1.915 


0.293 


0.6 


-0.963 


0.070 


4.281 


(-0.45,-0.34) 


0.364 


0.681 


1.912 


0.277 


0.7 


-0.970 


0.070 


4.298 


(-0.45,-0.34) 


-0.427 


0.682 


1.915 


0.278 


0.8 


-0.974 


0.066 


4.529 


(-0.45,-0.35) 


-1.097 


0.686 


1.919 


0.286 


0.9 


-0.976 


0.061 


4.885 


(-0.46,-0.35) 


-1.737 


0.693 


1.920 


0.297 


1.0 


-0.978 


0.056 


5.343 


(-0.46,-0.35) 


-2.390 


0.701 


1.923 


0.311 


1.1 


-0.979 


0.051 


5.879 


(-0.47,-0.35) 


-3.073 


0.708 


1.923 


0.324 


1.2 


-0.979 


0.046 


6.482 


(-0.48,-0.35) 


-3.796 


0.716 


1.916 


0.336 


1.5 


-0.981 


0.034 


8.768 


(-0.49,-0.34) 


-6.309 


0.738 


1.903 


0.375 


1.7 


-0.982 


0.028 


10.673 


(-0.50,-0.32) 


-8.317 


0.751 


1.887 


0.398 



13 For z eq = 1 and Q m = 0.3, any j3 in the interval (—1, —0.65) satisfies > (or rj^ > 0) and a > 0. However, the condition 

|<;| < 0.1 narrows down this interval to (—1, —0.978). 
14 See footnote 9 
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The second column in Table 1 gives the range of /3, for different values of z eq , that satisfy the requirements 
fcj^-j > 0, G > and \g\ < 0.1. This is illustrated in columns 3 and 4 by the positive values of i]k and (a/p m ), 
respectively. 

For these values, it turns out that the acceleration of the universe today is not an additional constraint or 
assumption. On the contrary it is a consequence of the model and the observational requirement |<?| < 0.1. The 
significant point here is that the results, which are presented in column 5, are consistent with current observations 
l)54[l. Besides, our model narrows down the possible values of q today 161(1 . 

For larger values of z eq , say z eq > 1.8, the requirements l|46() and (|49jl are automatically satisfied if (3 is selected 
in the range (—1,-0.983). However, the cosmic acceleration <|54[1 is not a consequence of the model (as it is for 
z e q < 1-8) but an observational constraint which restricts the values of (3 even further. 

Column 6 presents the dimensionless parameter rj ao . It changes sign for z eq 0.64. Therefore, there are two 
kinds of models: 

(i) Ever- expanding models. These are the ones with z eq < 0.64, for which r\ aa > 0. After a long matter-dominated 
phase of deceleration, the universe is accelerating today and will continue that way in the future 15 . 

(ii) Recollapsing models. For z eq > 0.64, we have i] ao < 0. After a long matter-dominated phase of deceleration, 
the universe is accelerating today, but deceleration comes back at some point in the future 16 . From l|64|l it follows 
that the expansion must come to an end before the universe starts contracting. Setting H = 0we can get a rec , the 
size of the universe at the moment of recollapse. Namely, 

- \ 3 , - x 3(/3+l) 

+(F-1) p- +2^=0. (68) 

Clearly, in the case under consideration the second term amply dominates over the first one. Therefore, from (|68|) 
we get 



F _ 1 \ 1/308+1) 



2<7 | 



0^-1. (69) 



The reason for the recollapse is that the four dimensional cosmological constant changes its sign sometime in the far 
future. Using (O and 1(23, we find 

a A(4) =o « W W+1) a rec . (70) 
As an example, we use the data for z eq = 0.7. We find 

a rec w 4.44 x 10 16 a, a A(4) w 0.71a rec . (71) 

The good news is that we are nowhere near the recollapse! 

The age of the universe t, in terms of the current value of the Hubble "constant" H, is presented in the seventh 
column. The main feature here is that the universe is older than the FRW counterpart for which Hi — 2/3. To put 
the discussion in perspective, we notice that in brancworld models (with A/4) = 0) the universe is younger than in 
the standard FRW cosmologies 17 . The age of the universe increases with z eq . As an example, if we take z eq = 100 
(for this value the crossover point is zt = 0.67), we find Hi = 0.96. 

Finally, the redshift of transition is specified in column 9. For the values considered here the crossover takes place 
at zt ~ 0.28 — 0.4. However, it can occur even earlier if the vacuum contribution starts to dominate for larger values 
of z eq , in which case the range for zt extends to 0.28 < zt < 0.68. This redshift interval coincides with the time of 
explosions of a number of SNe la known today. For example, SN 1988 at z = 0.31, SN 1992bi at z = 0.46, SN 1995K 

15 In order to avoid misunderstanding, let us mention that the equation ao(z e q, j3) = does have a solution for z eq = (0.4, 0.5, 0.6, 0.64) 
and p = (—0.67, —0.72, —0.83, —0.96), respectively. But these values of f) do not satisfy the physical conditions 1461 . 1491 and 1541 . 
Therefore they are excluded here. More precisely, the only models with erg = that satisfy physical conditions are those with z eq in the 
interval (0.642,0.643) and -1 < /3 < -0.967. 

16 For 0.64 < z eq < oo, using 1511 we find that equation ?(zt) = has two set of solutions. One for positive 2, which is given by l60l . 
and another one for negative z, namely —1 < zj> < —0.22. The second solution corresponds to the future, when the size of the universe 
is a > 1.28a. 

17 For example, for a dust-filled universe without cosmological term, we find ibrane = 0.363-ff — 1 in brane models, while tpjiyy = 2H~ 1 /3 
in FRW models. 



13 



at z = 0.48 and SN 1995 at z = 0.66. Therefore, they can provide crucial information to reduce the uncertainty in 
the transition between deceleration and acceleration. 

For completeness, we mention that the cosmological constant in the bulk, using (|26[1 and l|63[) . can be written as 

A (5) = -<^Hp m . (72) 

Thus, we can "predict" the value of A( 5 ) by means of measurements in AD. 



6 Discussion and conclusions 

An important distinction between general relativity and braneworld theory is that the cosmological term in the first 
is put by "hand" while in the second it is determined by the solution in the bulk through Israel's boundary conditions. 

A variable vacuum energy can generally occur as a consequence of embedding our universe as a brane in a five- 
dimensional bulk with non-static extra dimensions. There are a number of known solutions in 5D for which the 
vacuum density on the brane decreases as an inverse power of the scale factor, similar (but at different rate) to the 
power law in FRW-universes of general relativity. 

In this paper, we devoted our attention to spatially flat, homogeneous and isotropic, brane-universes where the 
vacuum density decays as in Q22[l, The model contains two parameters, viz., (3 and z eq . 

If (3 = 7, then A( 4 ) = 0, G is a universal constant and q — 2 — 3f2 m /2, in the dust dominated era. Thus fixing 
fl rn today also fixes q. Notice that q is positive for any physical value of f2 m , meaning that a brane- universe with 
constant vacuum energy must be slowing down its expansion. 

However, for f3 ^ 7, this is no longer so; the vacuum energy is now a dynamical quantity which changes this 
picture. In fact, G as well as A( 4 ) become functions of time and the deceleration parameter decreases from q 7 « 37 + 2 
at the beginning of the universe to q — > — 1 + 30 m /2 at the present time. Thus, for any matter with 7 > —2/3, q is 
positive at the beginning and negative today, because currently £l m < 2/3. 

Our model predicts that the transition from deceleration to acceleration occurred only recently, for zt < 0.68, 
but not later than zt ~ 0.28, regardless of the specific value of z eq . Therefore, early structure formation, from small 
density inhomogeneities, is not affected. 

If the domination of the vacuum is recent, i.e., z eq < 1.8, then the observed accelerated cosmic expansion is, not 
a condition but, a consequence of our model. The truth is that we have no observational or theoretical reasons to 
believe that vacuum started to dominate before z ~ 1. Consequently, we can safely declare that our model predicts 
the present acceleration of the universe. 

What is important here is that the predicted values for q are consistent with observations and allow us to narrow 
down the experimental uncertainty in the current data; from (15411 to 



q = -0.41 ± 0.09. (73) 

This may help in observations for an experimental/observational test of the model. 
Let us notice that an alternate way to write the effective density (|39l) is 

Peff=Pm + (74) 
Prn 

where n = (/3 + 1). This peculiar notation allows us to make contact between our ever-expanding models (those 
with do > 0) and Cardassian models. Indeed, for oo ^ Peff/% the second term in <|27|) can be neglected 18 . In this 
approximation, the generalized Friedmann equation (|27|) can be written as 

H 2 » A Pm + B P n m , (75) 

18 This is clearly satisfied for z eq as 0.4 during the recent epoch of accelerated expansion which started at 2 ^ 0.38. See Table 1 
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with n = (/3 + 1) and 

A = H 2 mV<? B = H 2 Vkr]* (F- 1) 

Pm ' P^ 

Equation l|75|) is similar to the one used in the so-called Cardassian models It differs from the usual Friedmann 
equation of general relativity by the addition of the extra term pj^. Therefore, based on (|75|) . one can interpret 
Cardassian expansion as the low energy limit of the brane model discussed here. Freese and Lewis £23 and Gondolo 
and Freese |41] suggest that the extra term may arise as a consequence of embedding our universe as a brane in extra 
dimensions. This is exactly what we have here. 

Nonetheless, there is a big difference in the understanding of the extra term. In the Cardassian model, by 
assumption, there is no vacuum contribution and the new term may come from some (yet unknown) modified 
Einstein equations. In our model the extra term is a manifestation of the variation of vacuum, in response to the 
time evolution of the extra dimension 

This difference leads to distinct requirements on n (or [3). In Cardassian models the most stringent requirements 
on n demand n < 0.4 (equivalently, f3 < —0.6). While, in our model the parameter f3 depends on Q m and is severely 
restricted by physical conditions fe?U > 0, G > and |<?| < 0.1. Indeed, for f2 OT = 0.3 we found —1 < (3 < —0.84 
(equivalently, < n < 0.16). 

In order to get another interpretation we write l|75|l as 



H 2 re H 2 



Q M (l + zf + n x (l + z) 3 ^ +1 A , (77) 



with wx — (3, CIm — VkVo-o an d — VkVcroi^ ~ !)■ I n this approximation F ps (??fcf7cr ) _1 , so that 17m + ~ 1- 
The above expression is simikar to the Friedmann equation in cosmological models with quintessence, where the dark 
energy component is characterized by the equation of state wx = Px/px- 

Thus, if the vacuum energy only started to dominate recently, at z eq ~ 0.4 — 0.6, then at low energy limit there 
is a correspondence between our ever-expanding as well as the Cardassian and quintessence models, although they 
are entirely different in philosophy. 

However, if the vacuum energy started to dominate earlier, say at z eq > 0.64, then do is negative and the quadratic 
term in (|27|l cannot be disregarded. As a consequence our model (for z eq > 0.64) is entirely different from Cardassian 
and quintessence models. In particular, it predicts that our universe will collapse in the future. 

We would like to finish this work with the following comments. The whole analysis in this paper is independent 
of any particular solution used in the five-dimensional bulk. This is a great virtue of brane-world models as noted 
at the end of Section 2. However, one could still ask whether the universe described here can be embedded in a five- 
dimensional bulk. The answer to this question is positive. For example, they can be embedded in five-dimensional 
"wave-like" cosmologies of the type discussed in . If in equation (38) of jHJ we take variable a as here in ill'L'l) , 
then the scale factor a for such wave-like models is governed by an equation which is identical to 1)25(1 in this paper. 

It is important to mention that the ratio ($/$) appears in different contexts, notably in expressions concerning 
the variation of rest mass ISHli electric charge [H3] and variation of the gravitational "constant" G |65 j .|66 j . 
Therefore, we have a scenario where the observed cosmic acceleration is just one piece in the dynamical evolution of 
an universe where the so-called fundamental "constants" are evolving in time. Braneworld scenario may provide us 
a theoretical framework to unify all these, apparently, separated phenomena as different consequences in 4D of the 
time evolution of the extra dimension. This is a new step toward understanding how the universe works. 

Note added in proof: While different alternative explanations are given to explain the current acceleration of 
the universe |21 | -|44 | . the dimming of the supernovae could be explained on the basis of axion physics |67j. 
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